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Numerical  Simulation  of  Quarry  Blast  Sources 


Summary _ 

This  tinal  report  summarizes  the  principal  results  of  numerical  investigations  of  the  quarry  blast  as  a 
seismic  source.  The  goal  was  to  derive  simple  but  physically  meaningful  models  that  could  be  used  to 
predict  regional  and  teleseismic  waveforms  for  large  industrial  blasts.  The  source  models  may  then  be 
used  to  develop  seismic  discrimination  strategies  for  the  identification  of  industrial  blasts  and  other 
natural  and  man  made  seismic  sources. 

Physical  models  are  presented  for  the  excitation  of  seismic  waves  from  quarry  blasts.  The  effects  of 
sp^l,  ripple  firing,  and  the  quarry  bench  on  the  excitation  of  seismic  waves  are  each  investigated. 
Theoretical  derivations  are  presented  for  physical  models  which  consist  of  forces  applied  to  the  free 
surface  to  model  spall  and  modifications  to  the  explosion  moment  tensor  to  model  the  effects  of  the 
non-planar  free  surface  of  the  quarry  bench.  Numerical  investigations  using  both  2D  and  3D  finite 
difference  calculations  are  used  to  verify  that  the  simplifying  assumptions  are  valid  for  the  band- 
widths  of  interest  It  is  concluded  that  the  spall  source  should  be  an  important  part  of  the  typical 
quarry  blast  and  that  horizontal  movement  of  material  away  from  the  quarry  face  should  provide  a 
source  of  additional  S  waves.  A  simple  numerical  model  for  the  effects  of  the  non-planar  quarry  face 
is  to  replace  the  non-deviatoric  explosion  moment  tensor  with  a  deviatoric  moment  tensor.  Therefore 
the  explosion  behind  a  quarry  face  will  appear  non-isotropic  and  radiate  more  S  waves  than  the  explo¬ 
sion  at  the  same  depth  in  a  layered  half-space. 

The  scenario  of  hiding  a  decoupled  or  partially  decoupled  underground  nuclear  explosion  in  the 
quarry  blast  is  considered  given  the  model  that  was  developed  for  the  quarry  blast.  It  is  concluded  that 
if  the  large  industrial  chemical  explosion  is  ripple  fired,  the  amplitude  of  the  seismic  signal  is  only 
10%  of  the  amplitude  of  a  blast  in  which  all  explosives  are  detonated  simultaneously.  Consequently, 
the  explosion  will  dominate  a  simultaneous  bomb  and  ripple-fired  quarry  blast  if  the  (tamped)  explo¬ 
sion  yield  exceeds  10%  of  the  total  quarry  yield.  Finally,  synthetic  regional  Pg  and  Lg  excited  by  the 
quarry  blast  model  are  compared  with  models  for  an  overburied  explosion,  a  normally  buried  explo¬ 
sion  with  spall,  and  a  shallow  crustal  earthquake.  Based  on  these  models  there  is  theoretical  justifica¬ 
tion  to  believe  that  frequency  dependence  of  Pg/Lg  ratios  and  Lg  spectra  may  contain  discrimination 
information. 

1.0  A  physical  model _ 

The  typical  large  quarry  blast  is  designed  to  frac¬ 
ture  and  move  large  quantities  of  rock  in  a  safe 
and  economical  manner.  To  accomplish  this, 
explosives  are  placed  at  scaled  burins, 
between  0.5  and  1  m  per  cube-root  Kg  of  explo¬ 
sives  (between  50  and  100  m/Kt^^^)  and  deto¬ 
nated  in  many  separate  holes.  Detonations  are 
timed  so  as  to  allow  the  rock  from  one  row  of 
charges  to  move  out  of  the  way  before  the  next 
row  of  charges  are  detonated.  In  this  way,  frac¬ 
ture  from  spallation  is  maximized.  The  blasting 
literature  (Langefors  and  Kihlstrom,  1963)  con¬ 
firms  that  the  spall  momentum  from  charges  det¬ 
onated  in  this  range  of  scaled  burden  is 
consistent  with  the  spall  momentum  predicted 
by  non-linear  finite  difference  simulations  of 
nuclear  explosives  in  the  same  range  of  scaled 
depth  of  burial,  4  *  10^®  Nt-sec  per  Kt  explosive 


Figure  l.Ol  A  typical  quarry  blast  consists  of 
multiple  detonations  of  up  to  several  tons  each 
behind  a  quarry  face. 


or  40  Nt-sec  per  Kg  of  explosive,  (see  Barker,  et 
al.,  1993a). 

See  Appendix  A  for  a  more  detailed  derivation  of 
the  equivalent  forces  required  to  model  the  spall 
from  a  quarry  face.  We  present  here  the  principal 
conclusions  of  Barker  et  al.  {1993a). 


Figure  1.02  A  physical  model  for  the  spall  asso¬ 
ciated  with  a  quarry  blast.  The  explosions  spall 
material  from  the  face  that  then  falls  to  the  quarry 
floor.  The  spalled  mass,  M,  may  fall  a  net  dis¬ 
tance,  h. 


To  describe  the  spall  source  for  a  quarry  blast, 
we  have  adapted  the  spall  model  of  Barker  and 
Day  (1990)  to  include  the  effects  of  horizontal 
throw  and  the  net  change  in  elevation  of  the  spal¬ 
led  mass.  Both  of  these  effects  are  found  to  make 
significant  contributions  to  the  seismic  source. 

The  equivalent  vertical  spall  force  is  given  by 

F,=  M[V^8(t)  +  (gtd  -V^)8(t-td)-g(H(t)-H(t-td))l  (EQ  1) 

and  the  horizontal  force  is  given  by 

Fx  =  M  Vxo[8(t)  -  5(t  -  td)],  (EQ  2) 

where  M  is  the  spall  mass,  and  W^q  are  the 
initial  vertical  and  horizontal  spall  velocities,  g 
is  the  acceleration  of  gravity,  and  t^  is  the  spall 
dwell  time, 

td  =  [Vzo  +  (VxO^  +  2h  g)  /  g,  (EQ  3) 


and  h  is  the  net  change  in  elevation  of  the  spall 
mass.  The  total  spall  momentum  is  P  =  M(V^2 

+  Vxo^)*^2  These  forces  are  applied  to  the  free 
surface.  Note  that  if  the  initial  horizontal  veloc¬ 
ity  is  zero  and  the  net  vertical  elevation  is  zero, 
h=0,  this  model  reverts  to  that  described  for  ver¬ 
tical  spall  in  Barker  and  Day  (1990).  The  model 
both  conserves  momentum  and  accounts  for 
energy  derived  from  the  net  loss  of  potential 
gravitational  energy.  The  only  variables  are  the 
initial  momentum,  the  total  mass,  and  the  direc¬ 
tion  of  initial  takeoff. 

The  time  functions  for  these  two  force  compo¬ 
nents  ate  plotted  in  Figure  1.03.  Obviously,  the 
takeoff  and  slap-down  phases  are  not  as  simple 
as  (telta  functions,  but  the  model  contains  the 
essential  elements.  The  time  functions  in  equa¬ 
tions  1  and  2  can  be  made  more  detailed  to 
reflect  more  complete  knowledge  of  the  spall 
ballistics. 
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Figure  1.03  The  dwell  time  is  determined  from 
the  magnitude  and  direction  of  the  takeoff  veloc¬ 
ity.  If  the  spall  mass  falls  a  net  distance,  h,  then 
there  is  a  contribution  from  the  gravitational  po- 

tential. 

We  have  computed  Green’s  functions  for  point 
forces  and  moment  tensor  sources  and  convolved 
them  with  the  spall  and  explosion  source  func¬ 
tions  expected  from  a  20  ton  explosion  buried  at 
a  depth  of  25  m.  The  results  are  shown  in  Fig¬ 
ures  1.04  and  1.05.  The  spall  is  a  significant 
source  compared  to  the  explosion. 


2.0  The  hide  in  quarry  blast  scenario 


jFigure  2.01  The  hide  in  quarry  blast  scenario. 

In  the  hide  in  quarry  blast  scenario,  a  large 
chemical  shot  designed  to  look  like  a  quarry 
blast  may  be  used  to  hide  an  overburied  nuclear 
explosion.  The  nuclear  explosion  may  be  par¬ 
tially  or  fully  decoupled.  We  have  examined  this 
scenario  under  the  assumption  that  the  shots  are 
fired  row-by-row  as  well  as  simultaneously.  In 
order  to  detonate  1  Kt  of  chemical  explosive 
using  standard  practices  one  scenario  would  be 
500  charges  of  2  metric  tons  each  fired  in  rows 
of  25  each.  Regional  synthetic  seismograms 
were  produced  at  a  distance  of  1000  km  for  a 
quarry  blast  of  20  by  25  2-ton  charges  and  for  an 
overburied  nuclear  explosion  of  1  Kt  at  a  depth 
of  300m.  This  scenario  is  one  of  many  such  pos¬ 
sible  configurations,  but  it  serves  to  illustrate  our 
primary  conclusions  and  is  consi.stent  with 
actual  bia.sting  practice. 

The  quarry  blast  synthetics  were  computed  using 
the  explosion  plus  spall  model  described  in  the 
previous  section.  Synthetics  were  examined  for 
several  cases,  including  pure  vertical  spall  and 
spall  with  an  initial  angle  of  45  degrees  to  the 
horizontal.  Green's  functions  were  computed 
from  0  to  5  Hz  using  a  reflectivity  code  and  con¬ 
volved  with  the  appropriate  source  time  func¬ 
tions.  The  spall  component  of  the  quarry  blast 
was  found  to  be  an  important  contribution  to  the 
quarry  blast  synthetic.  Example  synthetics  are 
shown  in  Figure  2.02  (see  also  Figure  2.03)  for 
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Figure  2.02  Vertical  component  syn¬ 
thetic  seismograms  (0.75-5  Hz)  at 
1(XX)  km  for  the  1  Kt  overburied 
bomb  and  the  1  Kt  ripple  tired  quarry 
blast.  The  sum  of  the  two  is  shown  at 
the  top.  Note  that  the  bomb  domi¬ 
nates  the  ripple  fired  quarry  blast. 
Note  the  different  Pg/Lg  ratio  of  the 
bomb  and  quarry  blast.  Note  that  the 
ripple  fired  quarry  blast  is  not  5(X) 
times  larger  than  the  single  2  ton 
quarry  blast  shown  at  the  bottom. 
This  is  because  of  the  extended  dura¬ 
tion  of  the  ripple  fired  quarry  blast. 
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Figure  2.03  Predicted  smoothed  Pg  spectra 
(top)  at  a  distance  of  1000  km  for  thel  Kt 
bomb  (dashed  line)  and  1  Kt  ripple  fired 
quarry  blast  (solid  line).  A  spectral  ratio  is 
shown  below  to  illustrate  the  enhanced  high 
frequencies  in  the  overburied  bomb. 
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Figure  2.0,4  From  top  to  bottom:  syn¬ 
thetic  seismograms  (0.75-5  Hz)  for  the 
I  Kt  overburied  explosion  (Bomb)  plus 
the  1  Kt  simultaneous  quarry  blast,  the 
1  Kt  overburied  explosion  only,  the  1  Kt 
simultaneous  quarry  blast  only,  and  a 
single  20  ton  shot.  Note  that  simulta¬ 
neous  detonation  of  500  2-ton  shots 
could  obscure  the  the  1  Kt  overburied 
explosion.  For  simultaneous  detona¬ 
tion,  the  1  Kt  quarry  blast  is  almost  ex¬ 
actly  500  times  larger  than  the  single  20 
ton  shot. 
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Figure  2.0.5  From  top  to  bottom:  synthet¬ 
ic  seismogram  spectra  (0-5  Hz)  for  the  1 
Kt  overburied  explosion  (Bomb)  plus  the 
1  Kt  simultaneous  quarry  blast,  the  1  Kt 
overburied  explosion  only,  the  1  Kt  si¬ 
multaneous  quarry  blast  only,  and  a  single 
20  ton  shot  Note  that  simultaneous  deto¬ 
nation  of  20x25  2-ton  shots  could  obscure 
the  the  1  Kt  overburied  explosion.  For  si¬ 
multaneous  detonation,  the  1  Kt  quarry 
blast  is  almost  exactly  500  times  larger 
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tte  0.75-5  Hz  bandwidth.  Because  the  ripple 
fired  quarry  blast  is  detonated  over  an  extended 
duration,  the  500  2  ton  blasts  generate  a  signal 
an  order  of  magnitude  smaller  than  an  amplitude 
500  times  the  amplitude  of  each  individual  2  ton 
blast.  In  order  to  hide  the  explosion,  the  500  sep¬ 
arate  charges  must  be  detonated  nearly  simulta¬ 
neously.  This  is  in  general  agreement  with  the 
results  of  Smith  ( 1 992). 

We  conclude  that  in  order  to  hide  an  overburied 
nuclear  explosion  (bomb)  with  a  seismic  yield  of 
1  Kt,  the  1  Kt  of  chemical  explosives  would  have 
to  be  detonated  nearly  simultaneously.  Even  so, 
the  overburied  explosion  may  still  be  visible  at 
high  frequencies  as  anomalous  high  frequency 
energy  in  the  quarry  blast  plus  bomb  seismo¬ 
gram  (see  Figures  2.04  and  2.05). 


3.0  The  effects  of  the  quarry  face _ 

3.1  A  theoretical  model 

A  theoretical  model  for  the  effects  of  the  quarry 
face  has  been  proposed  (Barker  et  al.  1993b). 

For  a  more  detailed  theoretical  discussion,  the 
reader  is  referred  to  Appendix  B.  The  essence  of 
this  model  is  that  for  a  moment  tensor  source 
located  behind  a  non-horizontal  surface  such  as 
a  quarry  face,  the  apparent  horizontal  couple 
perpendicular  to  the  face  is  effectively  reduced. 
Therefore  an  explosion  source  behind  the  quarry 
face  no  longer  looks  i.sotropic  when  compared  to 
an  explosion  in  a  half-space.  If  we  let  the  quarry 
face  be  aligned  along  the  Y  axis,  then  we  can 
summarize  this  model  by  saying  that  the  appar¬ 
ent  moment  tensor  component,  is  reduced 
from  that  of  the  actual  moment  tensor.  The 
apparent  moment  tensor  components  for  an 
explosion  behind  the  face  will  obey  the  relation. 

Mxx  <  Myy.  Myy=  ^zz- 

The  far-field  seismic  waveforms  are  then  com¬ 
puted  using  the  apparent  moment  tensor  compo¬ 
nents  and  plane-layered  Green's  functions  for  a 


model  without  the  quarry  face  .This  model  is 
valid  for  wavelengths  long  compared  to  the 
height  of  the  bench. In  order  to  illustrate  this 
result,  we  have  computed  several  numerical  sim¬ 
ulations  using  2D  and  3D  linear  finite  difference 
codes. 

3.2  2D  finite  difference  modeling 

2D  linear  finite  difference  calculations  have  been 
performed  using  reciprocity  to  examine  the  radi¬ 
ation  efficiency  of  a  point  explosive  source  in  the 
vicinity  of  a  quarry  face.  A  plane  P  or  Rayleigh 
wave  is  incident  upon  a  finite  difference  grid 
(Figure  3.2.1).  The  dilatation  is  then  monitored 
on  the  finite  difference  grid.  The  dilatation  for  a 
model  with  free-surface  topography  (like  a 
quarry)  is  then  compared  to  dilatation  in  a  half¬ 
space.  Likewise  displacements  monitored  in  the 
grid  are  used  to  infer  the  radiation  from  point 
sources.  The  result  is  a  map  of  the  efficiency  of 
radiation  by  sources  in  the  finite  difference 
model  relative  to  sources  in  a  half-space.  The 
results  are  contoured  in  Figure  3.2.2  for  a  plane 
P-wave  arriving  from  the  right. 


"igure  3.2. 1  Reciprocity  is  used  to  infer  the  far- 
field  radiation  of  P  and  Rayleigh  waves  by  mea¬ 
suring  dilatation  and  displacement  at  locations  in 
the  2D  grid  from  incident  plane  P  and  Rayleigh 
waves. 


The  results  shown  in  Figure  3.2.2  have  been 
found  to  be  insensitive  to  the  take-off  angle  of 
the  P-wave.  Similar  results  are  found  for  Ray¬ 
leigh  waves  arriving  from  either  the  right  or  the 


Source  Excitation  Near  a  Bench 


Figure  3.2.2  Excitation  effi¬ 
ciency  jf  a  point  .souice  as  a 
function  of  K«:ation  behind  a 
45  degree  quarry  face.  The 
moment  tensor  explosion 
.source  (top)  is  reduced  by  as 
much  as  5<)^  for  locations 
nehind  the  bench.  At  dis¬ 
tances  behind  the  bench 
greater  than  the  bench  height, 
the  explosion  source  relative 
to  the  half-space  is  only 
reduced  by  20%  or  less.  At  a 
location  just  below  the  toe  of 
the  quarry  face,  the  explosion 
source  is  enhanced  by  a  factor 
of  2  or  more.  The  vertical 
point  force  source  (bottom) 
is  hardly  affected  by  the  pres¬ 
ence  of  the  quarry  face.This 
analysis  is  valid  for  wave¬ 
lengths  much  longer  than  the 
bench  height. 


left.  In  all  cases,  the  excitation  from  a  point 
explosive  source  at  a  distance  behind  the  bench 
less  than  1/2  the  height  of  the  bench  is  reduced 
by  nearly  50%  relative  to  excitation  in  a  half¬ 
space. 

Simulations  using  incident  SV  waves  show  that 
SV  radiation  is  enhanced  by  the  presence  of  the 
bench  for  the  explosive  sources  relative  to 
sources  in  a  half-space.  The  explosion  behind 
the  bench  appears  more  deviatoric. 

The  2D  simulations  lead  to  several  predictions. 
First,  an  explosion  behind  a  quarry  face  will 
show  reduced  excitation  of  P-waves  and  Ray¬ 
leigh  waves  in  directions  perpendicular  to  the 
quarry  face  compared  to  an  explosion  in  a  half¬ 
space.  Point  forces  applied  to  the  quarry  face  and 
behind  the  quarry  face  will  not  be  so  affected. 
Half-space  Green’s  functions  may  be  used  to 
model  the  explosion  but  the  explosion  source 
should  be  replaced  with  a  non-isotropic  source 
with  Mxx  about  50%  of  Myy  =  M^x- 


A  corollary  to  the  model  that  M^x  does  not  equal 
Myy  is  that  the  explosion  behind  the  quarry  face 
will  excite  SH  waves  in  a  direction  45°  to  the 
bench  face  normal  and  there  will  be  an  azimuth- 
ally  dependent  radiation  pattern  to  the  P  and  SV 
waves.  Although  2D  modeling  can  not  explore 
this  conclusion  directly,  we  can  use  3D  finite  dif¬ 
ference  modeling  to  explore  the  azimuthal  .scat¬ 
tering  of  seismic  waves  generated  by  an 
explosion  behind  the  bench. 

3.3  3D  finite  difference  modeiing 

A  3D  finite  difference  calculation  was  per¬ 
formed  that  models  a  large  quarry  with  vertical 
faces.  A  point  explosive  source  was  located 
directly  behind  the  quarry  face.  The  geometry  of 
the  quarry  and  explosion  are  shown  in  Figure 
3.3.1.  The  Recursive  Grid  Refinement  (RGR) 
version  of  TRES-3D  was  used  to  perform  the 
calculation.  A  calculation  without  the  quarry  was 
computed  for  comparison. 


Figure  3.3.1  Geometry  of  the  3D  finite  difference 
calculation. 


Some  snap-shots  of  the  particle  velocity  wave- 
field  are  shown  in  Figures  3.3.2  and  3.3.3.  These 
snap-shots  illustrate  that  the  P  wave  amplitudes 
are  reduced  for  azimuths  perpendicular  to  the 
quarry  face.  These  results  are  consistent  with  the 
results  of  the  2D  simulations.  We  observe  non¬ 
isotropic  radiation  from  the  explosion  source 
behind  the  quarry. 


VarUeal  SacUon  Through  Tha  Sourea 
Exploaion  Sourea  Quarry 


Figure  3.3.2  Vertical  section  from  3D  finite  dif¬ 
ference  simulation.  The  P-wave  wavefront  is 
clearly  visible. 


The  3D  finite  difference  simulations  show  that 
the  radiated  P-wave  is  weakest  in  directions  per¬ 
pendicular  to  the  quarry  face,  as  if  the  M^x  com¬ 
ponent  of  the  explosive  source  were  reduced 
compared  to  the  Myy  component.  P-waves  radi¬ 
ated  parallel  to  the  quarry  face  are  not  signifi¬ 
cantly  different  from  the  half-space  case  (no 
quarry).  The  explosion  source  behind  the  quarry 
face  also  radiates  SH  waves. 

4.0  Discrimination 


An  important  problem  in  discrimination  is  the 
identification  of  quarry  blasts,  earthquakes,  nor¬ 
mally  buried  nuclear  explosions,  and  overburied 
or  decoupled  nuclear  explosions.  This  problem 
is  complicated  by  the  fact  that  future  events  of 
interest  will  occur  in  areas  where  we  have  little 
or  no  prior  experience  in  discriminating  explo¬ 
sions  and  earthquakes.  In  particular,  we  expect 
that  overburied  and  decoupled  explosions  will 
not  have  the  characteristic  spall  signature  of  nor¬ 
mally  buried  explosions.  Even  where  we  have 
experience  with  underground  nuclear  explo¬ 
sions  and  earthquakes,  we  have  somewhat  lim¬ 
ited  experience  with  overburied  explosions. 

Most  of  our  discrimination  experience  is  based 
on  normally  buried  explosions  compared  to 
earthquakes  or  on  large  quarry  blasts  compared 
to  earthquakes.  Therefore,  numerical  models  of 
seismic  sources  may  suggest  useful  discrimi¬ 
nants  for  new  areas  of  interest  provided  they  can 
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#  be  calibrated  to  regions  for  which  we  have  more 

extensive  discrimination  experience. 

In  Figure  4.01  we  show  comparisons  of  pre¬ 
dicted  Pg/Lg  spectral  ratios  for  four  different 
^  Mi^=4  seismic  sources  at  1000  km:  a  thrust 

earthquake  (depth  7.5  km),  a  ripple  fired  quarry 
blast,  a  normally  buried  explosion  with  spall, 
and  an  overburied  explosion  with  no  spall.  The 
Lg/Pg  ratio  is  .scaled  so  that  the  ratio  of  the 
^  earthquake  at  1  Hz  is  10  and  the  Lg  amplitude  at 

1  Hz  is  held  nearly  constant  (Ml=4).  The  details 
of  this  figure  depend  on  the  crustal  structures 
and  particularly  on  the  crustal  Q’s.  However  the 
qualitative  prediction  is  thrt  nt  frequencies  above 
0  2  Hz,  we  expect  that  the  earthquake  will  have 

enhanced  Lg/Pg  ratios  relative  to  those  at  1  Hz, 
the  normally  buried  bomb  with  spall  will  have 
the  lowest  Lg/Pg  ratios  at  3  Hz  relative  to  those 
at  1  Hz,  and  the  quarry  blast  and  overburied  or 
^  decoupled  bomb  will  be  intermediate.  These  pre¬ 

dictions  are  consistent  with  observations  of  Ben¬ 
nett,  et  al.  (1989).  In  studies  of  spectral  slopes, 
they  found 


that  the  ratio  of  Lg-to-P  fell  off  faster  with 
increasing  frequency  for  normally  buried  bombs 
than  earthquates  or  quarry  blasts.  In  Figure  4.02 
we  show  predicted  Lg  spectra  for  the  same  four 
sources.  The  spectra  are  all  normalized  to  1  at  1 
Hz.  The  quarry  blast  falls  off  faster  than  the 


earthquake  spectrum  also  consistent  with  obser¬ 
vations  of  Bennett  et.  al.  (1989).  Again  the 
details  of  this  plot  will  depend  on  the  particular 
crustal  model  chosen,  but  we  expect  that  some  of 
the  qualitative  results  will  apply  to  a  wide  range 
of  structures.  Because  of  the  spall  components  of 
the  quarry  blast  and  the  normally  buried  bomb 
we  expect  that  they  will  show  the  greate.st  spec¬ 
tral  decay  with  increasing  frequency.  These 
results  are  in  general  agreement  with  results  of 
Kim,  et  al.  (1993),  Wuster  (1993),  and  Hutchen- 
son  and  Herrmann  (1993),  for  quarry  blasts  and 
earthquakes. 


Figure  4.02  Normalized  Lg  spectra  for  four 
sources  at  1000km.  The  earthquake  is  the  most 
enriched  in  high  frequency  Lg. 


5.0  Conclusions 


We  have  proposed  models  for  the  “spall”  contri¬ 
bution  of  the  quarry  blast  and  the  effects  of  the 
quarry  face  on  the  explosive  source.  Detailed 
theoretical  discussions  are  contained  in  Appen¬ 
dices  A  and  B.  These  models  may  have  signifi¬ 
cant  implications  on  our  ability  to  identify  and 
discriminate  large  quarry  blasts  from  earth¬ 
quakes,  rock-bursts,  and  nuclear  explosions. 
These  models  for  radiation  from  the  spall  source 
and  explosion  source  behind  the  non-horizontal 
free-surface  could  be  tested  in  a  properly 
designed  experimental  program.  We  summarize 
our  specific  conclusions  below. 

•  Seismic  radiation  from  the  spall  part  of  a  typi¬ 
cal  quarry  blast  is  comparable  to  the  explo¬ 
sion  part  of  the  quarry  blast  source. 
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•  The  quarry  blast  source  is  a  band-limited  sig¬ 
nal  compared  to  the  overburied  bomb  due  to 
the  spall  signal  and  the  duration  of  the  distrib¬ 
uted  quarry  blast. 

•  Pg  spectral  values  rise  in  the  0.5-3  Hz  band 
for  the  overburied  bomb  but  decay  for  the 
quarry  blast  source  due  to  the  spall  contribu¬ 
tions  and  duration  of  the  quarry  blast. 

•  Quarry  blast  spectra  are  scalloped  primarily 
because  of  the  duration  of  the  source. 

•  The  explosion  will  dominate  a  simultaneous 
bomb  and  ripple  fired  quarry  blast  if  the 
(tamped)  explosion  yield  exceeds  10%  of  the 
total  quarry  yield. 

•  Near  simultaneous  detonation  of  the  entire 
quarry  blast  is  required  to  hide  an  overburied 
nuclear  explosion  of  comparable  yield. 

•  Seismic  amplitudes  are  insensitive  to  the 
direction  of  the  spall  as  long  as  there  is  mini¬ 
mal  elevation  change. 

•  A  net  elevation  change  comparable  to  the  bur¬ 
den  can  significantly  increase  the  seismic 
amplitudes  from  a  quarry  blast. 

•  2D  simulations  predict  that  P  and  Rayleigh 
radiation  from  a  point  explosion  behind  the 
quarry  face  may  be  reduced  by  a  factor  of50% 
due  to  the  free-surface  face. 

•  2D  simulations  predict  that  the  reduction  in 
seismic  signals  radiated  by  a  point  explosion 
behind  the  quarry  face  for  P-waves  radiated 
perpendicular  to  the  face  is  nearly  indepen¬ 
dent  of  slowness. 

•  2D  simulations  predict  that  seismic  signals 
from  the  point  force  equivalents  of  spall  are 
little  affected  by  the  free-surface  face. 

•  3D  simulations  predict  that  the  free  surface 
face  introduces  an  azimuthal  radiation  pattern 
and  generates  SH  waves. 

•  Source  models  for  earthquake,  quarry  blast 
with  spall,  overburied  explosion,  and  nor¬ 
mally  buried  explosion  with  spall  predict  that 
frequency  dependent  Lg  and  Lg/Pg  may  serve 
as  useful  discriminants. 

The  proposed  spall  model  for  seismic  radiation 

from  quarry  blasts  predicts  enhanced  SV  radia¬ 
tion  compared  to  the  point  explosion  source 

alone.  Furthermore,  the  directivity  of  spall  from 


a  non-horizontal  face  predicts  that  SH  waves 
will  be  generated.  Both  of  these  effects  make 
quarry  blasts  look  like  deviatoric  sources. 

Because  the  dwell  times  for  typical  quarry  blasts 
are  near  1  sec,  the  typical  quarry  blast  will  have 
a  peaked  spectra  and  it  is  possible  that  this  may 
lead  to  useful  regional  discriminants. 

The  recent  results  by  Kim,  etal.  (1993),  Hutch- 
enson  and  Herrmann  (1993),  and  Wuster  (1993) 
are  in  qualitative  agreement  with  our  source 
models;  small  earthquakes  are  enriched  in  high 
frequencies  relative  to  quarry  blasts. 
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Appendix  A 

Force  Equivalents  for  Spall 


where 


Consider  the  situation  that  a  volume  of  rock,  V, 
in  a  face  is  moved  along  a  surface,  £  (see  Figure 
A-1  below).  We  can  write  the  displacement, 

{/,  (X)  at  a  far-field  location,  X,  using  the  repre¬ 
sentation  theorem  (Aki  and  Richards,  1980), 

u.(X)  =  J<«:x,.(i^)6«^(ii)c.^^^(ii)G.^^(x.n) 

£ 

where  t/,  (X)  is  the  i’th  component  of  the  far- 
field  displacement  at  location  X,  (i^)  is  the 

normal  at  location,  n ,  on  the  surface  £,  Su^  (f))  is 
the  displacement  discontinuity  across  the  surface 

£,  C^*p,(^)  is  the  elastic  tensor,  and  G;p^(X,f|) 

is  the  q’th  derivative  of  the  i’th  component  of  the 
Green’s  tensor  for  displacement  at  X  given  a 
force  in  the  p’th  direction  at  location 


V  -  Spall  Vohiiiw 
X  -  Span  Surface 

Sv  -  Free  Surface  of  Spall  Volume 


and 

Aj.  =  jdl 
£ 

is  the  surface  area.  Su^  is  the  average  displace¬ 
ment  discontinuity  across  the  surface,  £.  Since 
G-^  satisfies  the  equation  of  motion  in  V,  we  can 
write, 

where  p  is  the  material  density  and  co  =  2nf  is 
the  angular  frequency.  If  we  now  integrate  over 
the  volume  V  and  convert  to  a  surface  integral, 

V  (I+Sv) 

where  (L  +  Sy)  is  the  combined  surfaces  enclos¬ 
ing  the  volume  V.  is  a  free  surface,  so  the  inte¬ 
gral  over  it  is  zero  and  we  have 

=  -(o'pVpG,* 

I  V 

Now  we  denote  by  G-^,  the  average  Green’s  ten¬ 
sor, 

»  V 


Note  that  from  reciprocity  we  can  interpret,  and 

r.,(^,X)  =  U;(Ti)C,.*„(ii).Gi^^(X,Ti) 


as  the  traction  in  the  k’th  direction  at  li  due  to  a 
displacement  in  the  i’th  direction  at  X. 

We  assume  that  Suj  (fi)  is  uncorrelated  with 
C.jp^(Tl)  ft)  on  the  surface,  £.  We  dis¬ 

cuss  this  assumption  further  below.  We  can  re¬ 
write  our  first  equation  as 

u.  (X)  =  (it)  (it)  G(p  ,  (t  ft) 


is  the  total  mass  in  the  volume,  V.  And  we  can 
finally  write, 

U.{X)  =  -<D  hUf^yG-^ 

where  the  seismic  source  is  a  force  proportional 
to  the  average  acceleration  time  history  of  the 
spall  mass.  Note  that  the  acceleration  is  averaged 
over  the  surface  of  the  spall  mass  while  the 
Green’s  function  is  averaged  over  the  volume  of 


the  spall  mass.  To  further  interpret  these  results, 
we  define, 

P^  - 

as  the  k’th  component  of  the  equivalent  spall 
momentum.  Then  the  equivalent  force  is  simply 
the  time  derivative  of  the  spall  momentum. 

It  is  easy  to  interpret  this  result.  As  the  material 
moves  away  from  the  face,  there  is  a  recoil  force. 
The  spall  temporarily  carries  away  with  it  the 
momentum  of  the  material  at  the  face.  The  mate¬ 
rial  then  falls  and/or  slides  down  the  face.  These 
resulting  impact  and  frictional  forces  bring  the 
material  to  a  halt.  The  recoil  and  impact  forces 
are  in  opposite  directions.  Unlike  the  explosion 
spall  from  a  contained  nuclear  explosion,  the 
forces  are  not  purely  vertical.  If  we  assume  the 
material  spalling  away  from  the  quarry  face 
behaves  as  a  ballistic,  then  we  can  write  the 
equivelent  forces  as  equations  1  through  3. 

We  now  briefly  discuss  the  conditions  that  per¬ 
mit  the  factoring  of  the  equation  to  arrive  at  our 
result.  This  assumption  that  is  uncorrelated 
with  on  the  surface  is  equivalent  to  the 

equation, 

Z 

where  the  overbar  indicates  the  average  over  the 
surface.  Note  that  from  the  reciprocity  theorem, 
we  can  interpret 

3G 

as  the  traction  in  the  k’th  direction  at  ^  due  to  a 

displacement  in  the  i’th  direction  at  X.  One  inter¬ 
pretation  is  that  the  displacement  discontinuity 
on  the  surface,  is  well  represented  by  its 

average,  6u^.  This  is  a  reasonable  assumption 
since  the  variations  in  the  displacement  disconti¬ 
nuity  will  excite  only  high  frequencies  waves 
and  cancel  out  in  the  far-fleld.  A  second  interpre¬ 
tation  is  that  this  equation  says  that  the  tractions 


on  the  surface  from  a  far-tield  source.  r,j  (n.  X) . 
are  well  represented  by  their  average  over  the 

surface,  X) .  Again  this  is  a  resonable 
approximation  for  wavelengths  longer  than  the 
size  of  the  quarry  face  that  has  spalled.  If  these 
conditions  are  met,  then  we  can  write  , 

fdI8«,\)C,  G  =  ^JdlC .  G 

J  i  J  Jip9  ‘P-9 

z  z 

which  leads  to  our  main  result. 


Appendix  B 

Effects  of  the  Quarry  Face 


We  are  interested  in  the  effective  radiation  pat¬ 
tern  from  an  explosion  source  behind  a  quarry 
bench.  The  geometry  is  shown  below.  The 
point  explosion  source  is  located  at  r?  i=  0.  a 
distance,  h,  to  the  left  of  the  vertical  bench,  h  is 
referred  to  as  the  burden  of  the  quarry  blast. 
Typical  quarry  blast  practice  is  to  choose  a  bur¬ 
den  between  30  percent  and  100  percent  the 
height  of  the  bench. 


b  •  8*nch  H«i0ht 

h  •  0urd«n  or  Sogreo  Oiftirw*  to  8«Aeh  Faeo 


where  oil,  022  and  033  are  stress  terms  at 

location  q  due  to  a  displacement  in  the  i'th 

direction  at  location  x  •  We  now  make  use  of 
this  reciprocal  problem,  where  Oij  is  the  stress 

tensor  at  due  to  a  displacement  at  Simi¬ 
larly,  we  will  write  uj  as  the  displacement  at  j) 
due  to  a  displacement  at  x  •  Dependence  on 
the  location  x  and  the  i'th  direction  of  motion 
is  suppressed  for  clarity  of  notation  in  the  fol¬ 
lowing  discussion. 

Now,  we  know  on  the  vertical  bench  (a  vertical 
free  surface)  that  <7jj  =0,  and  we  expand  the 

Til  dependence  of  Cu  in  a  Taylor's  expansion 
from  the  bench. 


The  seismic  displacement  at  a  location,  X , 
from  an  explosion  source  with  moment  tensor, 

Mjj  =  5ij  Mq  at  location  fj  may  be  written  as 


We  have  the  equation  of  motion  which  states 
that 


+^-|»7l=h=PUl(h) 


where  Uj(x)  i*  component  of  displace¬ 
ment  at  location,  x,K  =  A  +  2/i/3  is  the  bulk 
modulus  at  the  source  location,  Cijpq(n)  is  the 

elastic  tensor,  and  Gip.q(x,  rj)  is  the  q'th  deriv¬ 
ative  of  the  Green's  tensor  response  for  the  i'th 
component  of  motion  at  x  given  a  force  in  the 

p'th  direction  at  location  fj . 

Note  that  by  reciprocity, 

Ckkpq(n)Gip,q(X.  fj)  =  (TH(»i)+ <y22(^)+ <y33(^) 


together  with  the  free  surface  boundary  condi¬ 
tions  on  the  bench  (at  rji  =  b  ) 

hence  we  can  state  for  the  horizontal  equation 
of  motion  on  the  bench 

^l»7i=h=P“i('li=b) 

Therefore, 
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<^ll(»Ji=0)  =  -hpui(;7i=h) 


We  approximate  ui(r;i  =  h)  by  the  halfspace 
solution,  =  h),  and  we  can  write  . 

=0)  =  -hpu|(»S)(^,  =h) 

This  says  that  the  <5^  j  stress  from  a  displace¬ 
ment  source  at  location  x  is  proportional  to 
the  distance  from  the  bench,  h,  and  density,  p. 

How  much  is  the  explosion  behind  a  bench 
reduced  with  respect  to  an  explosion  in  a  half¬ 
space?  We  can  estimate  this  effect  by  writing 
the  expression  for  the  halfspace  Oj  j  compo¬ 
nent  of  stress  in  terms  of  the  halfspace  dis¬ 
placements. 


However,  it  is  clear  that  the  radiation  pattern 
from  an  explosion  behind  the  quarry  bench 
will  depend  upon  the  slowness  and  the  azimuth 
of  the  radiated  wave.  That  is  to  say  that 

r  =  r(k,h) .  A  complete  treatment  of  this  prob¬ 
lem  therefore  requires  the  solution  of  the  com¬ 
plete  scattering  problem  using  numerical 
techniques  as  in  the  section  3. 


<Tj((HS)  =(A  +2/i) 


+A 


Sr]\ 

- i - + - 

dr)2 


\ 


) 


To  estimate  the  ratio  ^  we  neglect  the 
second  term  to  get  . 

dill 

We  note  that  for  plane  P  waves  propagating  in 
the  r|  I  direction  we  can  make  the  approxima¬ 
tion  that 

dm 

so  we  have  that 

r=_£L!_  =  i*  =  kh 

~  pofoxi/HS)  “  « 

For.  h  =  10  m,  a  =  2000  m/s,  we  have  H  lOHz) 
=  0.3 1 ,  n.5Hz)  =  0. 1 5,  and  R I  Hz)  =  0.03 1 . 
This  would  indicate  that  the  P-waves  radiated 
perpendicular  to  the  bench  will  be  substan¬ 
tially  reduced  with  respect  to  the  half-space  for 
much  of  the  bandwidth  of  regional  seismo¬ 
grams. 


DISTWBUTION  LIST 

Michael  Browne 

NON-GOVERNMENT  CONTRACTORS 

Teledyne  Geotech 

3401  Shiloh  Road 

Garland,  TX  75041 

Prof.  Thomas  Ahrens 

Dr.  Lawrence  J.  Burdick 

Seismological  Lab,  252-2 1 

Woodward-Clyde  Consultants 

Div.  of  Geol.  &  Planetary  Sciences 

566  El  Dorado  Street 

California  Institute  of  Technology 

Pasadena,  CA  91125 

Pasadena,  CA  91 109-3245 

Dr.  Theodore  Cherry 

Dr.  T.  Bache,  Jr. 

Science  Horizons,  Inc. 

> 

Dr.  Thomas  J.  Sorena,  Jr. 

710  Encinitas  Blvd.,  Suite  200 

Science  Appl.  Int’I  Corp. 

10260  Campus  Point  Drive 

Encinitas,  CA  92024 

San  Diego,  CA  92 1 2 1 

Dr.  Kin  Yip  Chun 

Geophysics  Division 

Dr.  Peter  Basham 

Physics  Department 

Dr.  Robert  North 

University  of  Toronto 

Earth  Physics  Branch 

Geological  Survey  of  Canada 

Ontario,  Canada  M5S  1A7 

1  Observatory  Crescent 

Dr.  Paul  M.  Davis 

Ottawa,  Ontario,  Canada  KIA  0Y3 

Dept.  Earth  &  Space  Sciences 

> 

University  of  C^ifomia  (UCLA) 

Dr.  Douglas  R.  Baumgardt 

Dr.  Zoltan  Der 

Los  Angeles,  CA  90024 

ENSCO,  Inc. 

Prof.  Steven  M.  Day 

5400  Port  Royal  Road 

Department  of  Geological  Sciences 

\ 

Springfield,  VA  22151-2388 

Prof.  Jonathon  Berger 

San  Diego  State  University 

San  Diego,  CA  92182 

IGPP,A-025 

Dr.  Mark  D.  Fisk 

Scripps  Inst,  of  Oceanography 

Mission  Research  Corporation 

University  of  California,  San  Diego 

735  State  Street 

La  Jolla,  CA  92093 

Dr.  G.  A.  Bollinger 

P.O.  Drawer  7 1 9 

Santa  Barbara,  CA  93102 

Department  of  Geological  Sciences 

Prof.  Stanley  Flatte 

Viipnia  Polytechnic  Institute 

Applied  Sciences  Building 

1 

21044  Derring  Hall 

University  of  California 

Blacksburg,  VA  24061 

Santa  Cruz,  CA  9.5064 

The  Librarian 

Dr.  Roger  Fritzel 

Dr.  J.  Carter 

Pacific  Sierra  Re.search 

Dr.  Stephen  Bratt 

1401  Wilson  Blvd..  Suite  1 100 

Center  for  Seismic  Studies 

1.300  N.  17th  Street.  Suite  14.50 

Arlington.  VA  22209-3871 

Arlington.  VA  22209 

16 


Dr.  Holly  K.  Given 
Inst.  Geophys.  &  Planet.  Phys. 
Scripps  Inst.  Oceanography  (A-025) 
University  of  Califomia-San  Diego 
La  Jolla,  CA  92093 

Prof.  Hans-Peter  Baijos 
Institute  for  Geophysik 
Ruhr  University /Bochum 
P.O.  Box  102148 
4630  Bochum  1 ,  FRG 

Prof.  Donald  V.  Helmberger 
Seismological  Laboratory 
Div.  of  Geol.  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91 125 

Prof.  Eugene  Herrin 

Prof.  Brian  Stump 

Inst,  for  Study  of  Earth  and  Man 

Geophysical  Laboratory 

Southern  Methodist  University 

Dallas,  TX  75275 

Prof.  Bryan  Isacks 
Prof.  Muawia  Barazangi 
Cornell  University 
Department  of  Geological  Sciences 
SNEE  Hall 
Ithaca,  NY  14850 

Prof.  Lane  R.  Johnson 
Prof.  Thomas  V.  McEvilly 
Seismographic  Station 
University  of  California 
Berkeley,  CA  94720 

Dr.  R.  Kemerait 
ENSCO,  Inc. 

445  Pineda  Court 
Melbourne.  FL  32940 

Prof.  Brian  L.  N.  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
G.P.O.  Box  4 
Canberra  2601.  Australia 


Dr.  Richard  LaCoss 
MIT-Lincoln  L^ratory 
M-200B 
P.O.  box  73 

Uxington,  ?MA  02173-0073 

Prof.  Fred  K.  Lamb 
University  of  Illinois 
Dept,  of  Physics 
1110  West  Green  Street 
Urbana,  IL  61801 

Prof  Charles  A.  Langston 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Prof.  Thome  Lay 
Dr.  Susan  Schwartz 
Institute  of  Tectonics 
Earth  Science  Board 
Univ.  of  California,  Santa  Croz 
Santa  Cmz,  CA  95064 

Prof.  Arthur  Lemer-Lam 

Prof.  Paul  Richards 

Prof.  C.  B.  Scholz 

Lamont- Doherty  Geol.  Observatory 

Columbia  University 

Palisades,  NY  10964 

Dr.  Manfred  Hanger 

Fed.  In.st.  forGeosci.  &  Nat’I  Res. 

Postfach  510153 

D-.3000  Hanover  51.  FRG 

Dr.  Peter  Marshall 

Procurement  Executive 

Ministry  of  Defen.se 

Blacknest.  Brimpton 

Reading  FG7-4RS.  United  Kingdom 

Dr.  Randolph  Martin.  Ill 
New  England  Research.  Inc. 

76  Olcott  Drive 

White  River  Junction.  VT  05001 


17 


Dr.  Bernard  Maseinon 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  ftiris,  France 

Dr.  Gary  McCartor 
Prof.  Henry  L.  Gray 
Department  of  Physics 
Southern  Methodist  University 
Dallas,  TX  75275 


E)r.  Frode  Ringdal 
Dr.  Svain  Mykkeltveit 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  KJeller,  Norway 

Dr.  Wilmer  Rivers 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314 


Dr.  Keith  L.  McLaughlin  Dr.  Richard  Sailor 

S-CUBED,  A  Div.  of  Maxwell  Laboratories,  Inc.  TASC,  Inc. 

P.O.  Box  1620  55  Walkers  Brook  Drive 

La  Jolla,  CA  92038- 1 620  Reading,  MA  0 1 867 


Dr.  Pierre  Mecheler 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris,  France 

Prof.  Bernard  Minster 
Prof.  John  Orcutt 
IGPP,A-025 

Scripps  Inst,  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Prof.  Brian  J.  Mitchell 
Dr.  Robert  Herrmann 
Dept.  Earth  &  Atmos.  Sci. 

St.  Louis  University 
St.  Louis,  MO  63156 

Mr.  J.  Murphy 
Maxwell  Laboratories,  Inc. 
S-CUBED  Reston  Geophysics  Office 
1 1 800  Sunrise  Valley  Dr.,  Suite  1212 
Reston,  VA  2209 1 


Prof.  Charles  G.  Sammis 

Prof.  Kei  Aki 

Center  for  Earth  Sciences 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089-0741 

Prof.  David  G.  Simpson 
Lamont-Doherty  Geol.  Observatory 
Columbia  University 
Palisades,  NY  10964 

Dr.  Stewart  W.  Smith 
Geophysics  AK-50 
University  of  Washington 
Seattle,  WA  98195 

Prof.  Clifford  Thurber 
Prof.  Robert  P.  Meyer 
University  of  Wisconsin-Madison 
Dept,  of  Geology  &  Geophysics 
1215  West  Dayton  Street 
Madison,  WS  53706 


Dr.  Jay  J.  Pulli 
Radix  Systems,  Inc. 

2  Taft  Court,  Suite  203 
Rockville.  MD  208.50 


Prof.  M.  Nafi  Toksoz 
Prof.  Anton  Dainty 
Earth  Resources  Lab 
Massachu.setts  Inst,  of  Technology 
42  Carleton  Street 
Cambridge,  MA  02142 


18 


Prof.  Terry  C.  Wallace 
Dept,  of  Geosciences 
Building  77 
University  of  Arizona 
Tucson,  AZ  85721 

Dr.  William  Wortman 
Mission  Research  Corporation 
735  State  Street 
P.O.  Drawer  719 
Santa  Barbara,  CA  93102 


U,S.  GOVERNMENT  AGENCIES 


Mr.  Alfred  Lieberman 
ACDAA^l-OA,  Room  5726 
320  21st  Street,  N.W. 
Washington,  D.C.  2045 1 

Colonel  Jerry  J.  Perrizo 
AFOSR/NP,  Building  410 
Bolling  AFB 

Washington,  DC  20031-6448 

Dr.  R.  Blandford 
AFTAOCSS 

1300  N.  17th  Street,  Suite  1450 
Arlington,  VA  22209-3871 

AFTAaCA 

STINFO 

Patrick  AFB,  FL  32925-6001 

Dr.  Frank  F.  Pilotte 

HQ  AFTAC/TT 

Patrick  AFB.  FL  32925-6001 

Katie  Polay 
CIA-ACIS/TMC 
Room  4X16NHB 
Washington,  DC  20515 

Dr.  L.  Turnbull 
CIA-OSWR/NED 
Washington.  DC  20505 


Dr.  R.  Alewine  and  Dr.  A.  Ryall 
ARPA/NMRO 

3701  N.  Fairfax  Drive,  Suite  717 
Arlington,  VA  22203-1714 

ARPAyOASB/Librarian 
3701  Fairfax  Dr.,  Suite  717 
Arlington,  VA  22203-1714 

Dr.  Dale  Glover 
DIA/DT-IB 
Washington,  DC  20301 

Dr.  Michael  Shore 
Defense  Nuclear  Agency/SPSS 
6801  Telegraph  Road 
Alexandria,  VA  22310 

Mr.  Max  A.  Koontz  (DP-5) 

Department  of  Energy 
Forrestal  Building 
1000  Independence  Avenue,  S.W. 
Washington,  DC  20585 

Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  VA  22314 

Dr.  John  J.  Cipar,  PL/LWH 
Phillips  Lab/Geophysics  Directorate 
Hanscom  AFB,  MA  01731 

Mr.  J.  Lewkowicz,  PL/LWH 
Phillips  Lab/Geophysics  Directorate 
Hanscom  AFB,  M A  0 1 73 1 

Phillips  Lab  (PL/XO) 

Hanscom  AFB,  MA  01731 

Dr.  J.  Hannon  and  Dr.  K.  Nakanishi 
Lawrence  Livermore  Natn'l  Laboratory 
PO.  Box  808 
Livermore,  CA  94550 

Office  of  the  Secretary  of  Defense 
DDR&E 

Washington.  DC  20330 


19 


Eric  Chael 
Division  924 1 
Sandia  LabonUory 
Albuquerque,  NM  87175 

Dr.  William  Leith 
U.S.  Geological  Survey 
Mail  Stop  928 
Reston,  VA  22092 

Dr.  Robert  Masse 
Box  25046,  Mail  Stop  967 
Denver  Federal  Center 
Denver,  CO  80225 

Directed  Energy  Effects 
PL/WSM 

Kirtland,  NM  87117-5776 


